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Abstract. Light gravitino productions in association with a neutralino (selectron) in e + e~ (e~ 7 ) collisions 
are restudied in a scenario that the lightest supersymmetric particle is a gravitino and the produced 
neutralino (selectron) promptly decays into a photon (electron) and a gravitino. We explicitly give the 
helicity amplitudes for the production processes by using the effective goldstino interaction Lagrangian, 
and present the cross sections with different collision energies and mass spectra. We also examine selection 
efficiencies by kinematical cuts and beam polarizations for the signal and background processes, and show 
that the energy and angular distributions of the photon (electron) can explore the mass of the t-channel 
exchange particle as well as the mass of the decaying particle at a future e + e~ (e _ 7 ) collider. 



1 Introduction 

Gravitinos are spin-3/2 superpartners of gravitons in local 
supersymmetric extensions to the Standard Model (SM). 
Since the gravitino becomes massive via the super-Higgs 
mechanism, its mass is related to the scale of supersym- 
metry (SUSY) breaking as well as the Planck scale like 



%/2 



(M SU sy)7^pi- 



(1) 



This implies that the gravitino can take a wide range of 
mass, depending on the SUSY breaking scale, from eV 
up to scales beyond TeV, and provide rich phenomenol- 
ogy in particle physics as well as in cosmology pQ. While 
the interactions of the helicity ±3/2 components of the 
gravitino are suppressed by the Planck scale, those of the 
helicity ±1/2 components are suppressed by the SUSY 
breaking scale if the gravitino mass is much smaller than 
the energy scale of the interactions, due to the goldstino 
equivalence theorem, and can be important even for col- 
lider phenomenology. 

Gravitino productions in association with a SUSY par- 
ticle are known processes which become significant at col- 
liders when the mass of the gravitino is very light as 
m.3/2 ~ 0(1O~ 2 eV) or less, since the cross sections are 
inversely proportional to the square of the gravitino mass 



a oc l/m 2 .^ 



(2) 



Such a very light gravitino is suggested by the context 
of no-scale supergravity [3J|3] and some extra-dimensional 
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models [I], while typical gauge-mediated SUSY breaking 
(GMSB) scenarios expect a mass of 1 eV-10 keV pQ. Sev- 
eral studies on the associated gravitino productions have 
been performed so far, for instance, Xi~G productions in 
e + e~ [5H71I8] and hadronic [5] collisions, e-G productions 
in e 7 collisions [5], and g-G [TUEIJEEi] and q-G [TEEl] 
productions in hadronic collisions. When the associated 
SUSY particle is the next-to-lightest supersymmetric par- 
ticle (NLSP) and promptly decays into a SM particle and a 
LSP gravitino, the above production processes lead to par- 
ticular collider signatures, such as 7+$, e+$, and jet+j£, 
where the missing energy is carried away by two graviti- 
nos, and these signals set mass bounds on the gravitino 
and the other SUSY particles. The current experimental 
bound on the gravitino mass from the single-photon plus 
missing-energy signaQ in x?-G associated productions is 
given by the LEP experiment as a function of the neu- 
tralino and selectron masses [T3], e.g. 



m 3 / 2 



> 10 -5 eV 



(3) 



for m%p 



140 GeV and m s = 150 GeV. We note that 
the Tevatron also set a similar bound on the gravitino 
mass for the 7 + f [T3] and jet±j£ [TS] channels, where it 
is assumed, however, that all SUSY particles except the 
gravitino are too heavy to be produced on-shell [T6] . 

While the previous searches for gravitino productions 
have been somewhat restricted due to limitations of simu- 
lation tools, in the recent paper |17j implementation of the 



We note that a two-photon plus missing-energy signal, 
where the two photons come from two neutralino decays, does 
not provide any constraint on the gravitino mass. 
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spin-3/2 gravitino in MadGraph/MadEvent (MG/ME) [TW2TH 
|2T] was reported, where new HELAS (HELicity Amplitude 
Subroutines) [TB] codes were introduced to calculate he- 
licity amplitudes with massive spin-3/2 gravitinos. They 
are implemented in such a way that amplitudes with exter- 
nal gravitinos can be generated automatically by MG/ME0 
Since goldstinos appear as the longitudinal modes of mas- 
sive gravitinos and their interactions become dominant 
over the transverse modes in high-energy processes, two 
of the authors also implemented effective goldstino inter- 
actions [22] as an alternative to the gravitino code. 

In this paper, we revisit the following two processes by 
using the gravitino implemented MG/ME mentioned above. 
First, we study associated gravitino productions with a 
neutralino which promptly decays into a photon and a 
gravitino in e + e~ collisions, 

e+e- -> x°iG -> 7GG, 

in the context of a neutralino NLSP with a gravitino LSP. 
In order to investigate the production cross section and 
distributions of the photon in detail, we explicitly give 
the helicity amplitudes for the production process by using 
the effective goldstino interaction Lagrangian, and present 
the cross sections with different mass spectra and different 
energies especially for a future linear collider. We also ex- 
amine selection efficiencies by kincmatical cuts and beam 
polarizations for the signal and SM background processes, 
and show that the energy and angular distributions of the 
photon coming from the neutralino decay can explore the 
mass of the i-channel exchange selectrons as well as the 
mass of the decaying neutralino. 

Second, we consider gravitino productions in associa- 
tion with a selectron which subsequently decays into an 
electron and a gravitino at an ej collider, which is an op- 
tion at a future linear collider |23j . 

e~7 -t e~G -> e~GG, 

in a slepton co-NLSP scenario with a gravitino LSP. We 
present the explicit helicity amplitudes for the produc- 
tion process, and discuss the mono-electron plus missing- 
energy signal, including the Compton back-scattered pho- 
ton energy spectrum 24, 2S] for incident photons. While 
the heavy-mass limit for all SUSY particles except grav- 
itino and selectron are assumed in Ref. [9], we take into 
account the i-channel intermediate ncutralinos and show 
a possibility to determine their mass in the signal distri- 
butions. 

We note in passing that all the helicity amplitudes we 
present are easily applicable to qq — > gG and to qg qG 
subprocesses for hadron colliders. 

The paper is organized as follows: In Sect.[5]neutrahno- 
gravitino productions in electron-positron collisions are 
considered, and in Sect. [3] selectron-gravitino productions 
in electron-photon collisions are studied. Sect.|3]is devoted 
to our summary. In Appendix[X]we give the effective gold- 
stino interaction Lagrangian relevant to our study, and in 
Appendix [B] we briefly mention neutralino decays into a 
photon and a gravitino. 

2 The spin-3/2 functionality is available in MG/ME V4.5 [21] . 



2 Neutralino-gravitino production in e + e 
collisions 

In this section, we consider a scenario of a neutralino 
NLSP with a gravitino LSP, and study associated grav- 
itino productions with a neutralino which promptly de- 
cays into a photon and a gravitino in e + e~ collisions, 

e+e" x\G -> 7GG, (4) 

leading to a mono-photon plus missing-energy signal. 

2.1 Helicity amplitudes 

Here we present the helicity amplitudes explicitly for the 
production process: 

(5) 

where the four-momentum (p^ and helicity (Aj = ±1) of 
each particle are defined in the center-of-mass (CM) frame 



Diagrams by MadGraph e- e+ -> n1 gld 




graph 5 graph 6 



Fig. 1. Feynman diagrams for the Xi-G production in e + e~ 
collisions, generated by MadGraph [22] . Nl, gld, el, and er de- 
note a lightest neutralino, a gravitino, a left-handed selectron, 
and a right-handed selectron, respectively. 
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Table 1. The reduced helicity amplitudes M\,\ 3 \ l for e x 



of the e + e~ collisions. Throughout our study only the he- 
licity ±1/2 components of the gravitino, i.e. goldstinos, 
are considered. In the massless limit of e , one can find 
that all the amplitudes are zero when both the electron 
and the positron have the same helicity, or Ai = A2. In 
addition, for the Ai = +1 (Ai = — 1) case, only the right- 
handed (left-handed) selectron can contribute to the to- 
tal amplitudes. Therefore, the helicity amplitudes for the 
above process can be expressed as the sum of s-, t-, and 
u-channel amplitudes: 

M XMX4 = M{ MXi + M\ MXi + Ml MXi (6) 

with A = Ai = — A2, where each amplitude with A = +1 
(A = —1) corresponds to the Feynman graph 3+4, 2 (1), 
and 6 (5), respectively, in Fig. Q] 

We first present the amplitudes based on the effective 
goldstino interaction Lagrangian, given in Appendix [21 in 
the usual four-spinor basis: 



iM 



v(p 2 ,-\)Yu( Pl ,\) 



a,a 3 a 4 - 2 V6Mpim 3/2 s 

x u(p 3 ,\ 3 )ty 3 +p i ,-f ll ]v(p i , A 4 ), (7a) 
-VSeCf^ml 1 



V3Mpim 3 /2 t ^ m \ x 
x u(p 3 , A 3 )u(pi, A) v(p 2 , -X)v(pi, A 4 ), (7b) 
-yf2eCf^m\ 1 

x u(p 4 , A 4 )u(pi, A) v(p2, -X)v(p 3 , A 3 ), (7c) 



^A,A 3 A 4 



where M P1 = M Pi /\/8iT ~ 2.4 x 10 18 GeV is the reduced 
Planck mass, mg ± denotes the right-/left-handed selectron 
mass for notational convenience, and 



respectively^ The couplings related to the neutralino mix- 
ing denned by X l = U l3 x° in the X = (B,W 3 ,H°,H°) 
basis, where Uij is taken to be real, are 



CTC =Uu cos8 w + U 2 i sin 6w, 
C Zxi — — Uu sin 8w + Ui% cos 9w , 
U 2i , „ ? U u 



sin t)\y 



(10) 



with the SU{2) charge T| and the 17(1) charge Fj for 
e_)_/_(= e.R/1,). Here, for simplicity, we assume the lightest 

neutralino as a pure gaugino, which makes the G-H® U -Z 
couplings irrelevant to our study. It should be noted that 



and C^ Xl are related with each other as 



(11) 



for t/s 3> mz] this is always the case in the following 
discussions. 

To present the explicit helicity amplitudes, let us now 
define the kinematical variables of the process (J5J) in the 
e + e~ laboratory frame as 



Pi = 


2 


;i, 0,0,1), 


P2 = 


2 


;i, 0,0,-1), 


Ps = 




(1 + ^,0 sin 9,0,13 cos 9), 


n = 




(1 - ^,-/3 sin 6, 0,-/3 cos 0), (12) 


with = 


1- 


-m?o /s. Throughout our study we neglect the 



gravitino mass, except in the gravitino couplings. 

For notational convenience we define the reduced he- 
licity amplitudes, M., as 



C x = C 1Xl 



s — m z + imzTz 



9x 



c z xl 



(8) »A^a,a 3 a 4 = r-— e \/0sM XMX4 

v6Mpim 3 /2 



(13) 



with Z-boson couplings to right- and left-handed charged 
leptons, 



sin 9w — 1 + 2 sin 9w 

9+ = — and -9- = TT~- — n n — ' 

2 sin Uw cos 9w 



cos 9 W 



(9) 



and these are presented in Table [T] The following features 
of the amplitudes are worth noting: 



3 Strictly speaking, the Z-exchange amplitude in ()7a|) is valid 
only for y/s 3> mz since massless gauge bosons are assumed in 
the effective Lagrangian (|38|l . 
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1. As mentioned before, for the A = +1 (A = —1) case 
only e + (e_ ) can be exchanged in the t- and u-channel 
amplitudes, and all the amplitudes are zero for Ai = 
A 2 . 

2. The overall angular dependence is dictated by J = 1 
d functions as 



M x ,x 3 x i oc d{ (Xa _ Xi)/2 (6). 



(14) 



A4 S and J\A t,u interfere subtractively with each other; 
especially for the A3 = A4 case they almost cancel 
in the wide range of the parameter space, and hence 
the amplitudes with A3 = — A4 are dominant for the 



most of the cases except for the 



nrie ± region. 



We note that in the very high-energy region the am- 
plitudes with A3 = A4 become important since A4 S 
becomes dominant, making those amplitudes be pro- 
portional to y/s while the amplitudes with A3 = — A4 
are independent of \fs; in that region the cross section 
does not depend on the selectron masses but on the 
produced neutralino mass. 

4. For the A3 = — A4 case, in the threshold region, where 
t,u = —s0(l =F cos0)/2 — !• 0, an additional j3 can be 
extracted from the reduced amplitudes due to C± ~ 
-Cl* 1 in (HU). Together with p 1 ' 2 in {T3J), the am- 
plitudes are proportional to /3 3 / 2 . Therefore, including 
the phase space factor j3, the threshold excitation of 
the total cross section is given by [oTT UHpl 



a cx /3 4 . 



(15) 



5. A4 f and M. u depend on the selectron mass and be- 
come larger as the selectron mass increases, while A4 S 
is independent of mg. 

We note that our helicity-summed amplitude squared 
agrees with Eq. (28) in [8 for the photino case, and also 
with Eq. (3) of for the gluino associated process qq — > 
gG after substitutions for the masses and the couplings as 



-9sT a , 



le R/L 



(16) 



Moreover, our analytic amplitudes are checked numer- 
ically for each helicity combination by using the grav- 
itino/goldstino code in MG/ME [T7Il2"2] . 



2.2 Cross sections and kinematical distributions 

Let us now present the total cross sections and the kine- 
matical distributions for the production process (0. The 
initial-helicity (A) dependent cross section is given by 



^3,4 



(17) 
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Fig. 2. Total cross sections of associated gravitino produc- 
tions with a bino-/wino-like neutralino in e + e~ collisions, 
e + e~ — > XiG, for m 3 / 2 = 1CP 13 GeV as a function of the col- 
lision energy. The neutralino and selectron masses are fixed at 
300 GeV and 800 GeV, respectively. The initial-helicity depen- 
dent cross sections ax are shown by a dashed line for A = +1 
and a dotted line for A = — 1. The threshold region is enlarged 
and the hypothetical dependence a oc /3 and /3 s with the same 
coefficient is also shown. 



with the two-body phase space factor d<S>i- c un poi = (° + + 
<x_)/2 is the usual spin-summed and averaged cross sec- 
tion. 

Figure [5] shows total cross sections of the gravitino 
productions associated with a bino-/wino-like neutralino 
in e + e~ collisions as a function of the CM energy ■^/s, 
where the neutralino and selectron masses are fixed as 

800 GeV. It should be 



m^o — 300 GeV and mg. 



me 



stressed that the cross section scales with m 3 y 2 , and we fix 

~ 13 GeV in our study so that 
10 3 ) 



the gravitino mass m 3 / 2 = 10 

the production cross sections are around ©(lO 1 ~ 10 3 ) fb. 
In the figure the threshold region for the unpolarized bino- 
like neutralino cross section is enlarged, and one can see 
that the production cross section is strongly suppressed 
as shown in (|15[) . in contrast to the threshold excitation 
for the standard fermion (cx f3) and the scalar (oc (3 3 ) pair 
productions This is one of the particular signatures 
for the associated gravitino productions. 

For the case of the bino-like neutralino, or \Un\ ~ 1 
and \U2i \ ~ in (ITU1) . the cross section with right-handed 
electrons (0+) dominates the one with left-handed (<r_). 
For the heavy selectron case the t- and u-channel contri- 
butions are dominant, and therefore the ratio of the A- 
dependent cross sections is roughly given in terms of the 



Note that /3 is defined as (1 - m? /s) 1/2 in Refs. [7118] . 
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3. Total cross sections of e + e~ — > Xi@ a ^ \fs = 500 GeV 
1 TeV for m 3 / 2 = 10~ 13 GeV as a function of the neu- 
tralino mass. The selectron masses are fixed at 400 (solid), 800 
(dashed) and 1200 (dotted) GeV, respectively. 



Xi~e-e± couplings as 



2 



2(7 



unpol 



c; xi | 2 + \ct 



eXl 12 



(18) 



The bino case gives <7+/2<7 unp oi ~ 0.8, which one can ob- 
serve in Fig. [5] On the other hand, for the case of the 
wino-like neutralino, or |Z7ii| ~ and \U2i\ ~ 1, the right- 
handed cross section vanishes, i.e. cr unpo i = <r_/2. One can 
conclude that the Xi'G production process with polarized 
electron beam can explore the neutralino mixing. The de- 
tailed study for various neutralino mixing has been done 
in Rcfs. [7,8 , while we assume a bino- like neutralino in 
the following analyses for simplicity, which is often the 
case in GMSB; see, e.g. Fig. 4 in [27], 

In Fig. [3J the neutralino-mass dependence of the cross 
sections is shown for y/s — 500 GeV and 1 TeV. Due to the 
threshold behavior in (TT5j) , the cross sections are strongly 
suppressed as the neutralino mass is approaching the col- 
lider energy(f| It should be emphasized here that the cross 
section is quite sensitive to the mass of the t, w-channel 
intermediate selectrons, even if the collider energy cannot 
reach them [TJIH] ■ The heavier selectron exchange increases 
the cross section since the t, u-channel amplitudes are pro- 
portional to the selectron mass squared as one can see in 
Table [TJ We also note that, however, the goldstino cou- 
plings become too strong at some point for heavy selec- 
trons to perform the reliable perturbative calculations. 



5 For the case of mg ± = 400 GeV at yfs — 1 TeV, the cross 
section is not so strongly suppressed as /3 4 . This is because in 
this parameter region the contributions from the amplitudes 
with A3 = A4 are significant and these amplitudes do not pro- 
vide an additional suppression factor f3. 




Fig. 4. Normalized angular distributions of the neutralino in 
e+e" x°G at y/a = 500 GeV (left) and 1 TeV (right) for 
m^o = 300 GeV, where the selectron masses are taken to be 
400 1 , 800 and 1200 GeV. 



Before taking the neutralino decay into account, we 
discuss the angular distribution of the produced neutralino 
since the \i lG decay is isotropic (see Appendix iBl) and 
hence the photon distribution is given by purely kinemati- 
cal effects of the decaying neutralino. In Fig.^J the normal- 



ized cos 9 distributions of the neutralino in 



at v/i = 500 GeV (left) and 1 TeV (right) are shown for 
m^o = 300 GeV. One can find that not only the total cross 
section as shown in Fig.[3]but also the angular distribution 
is quite sensitive to the mass of the t, it-channel intermedi- 
ate selectrons [5] . When the selectron mass is close to the 
neutralino mass, the cross section is suppressed around 
I cos#| = 1 since a cancellation occurs between Ai s and 
M t,u for the A3 = — A 4 case due to (t — m|) = — s for 
cos (9 = —1 and (it — m|) = — s for cos 9 — 1; see also 
Table[TJ For the heavy selectron case, on the other hand, 
the neutralino tends to be produced to the forward and 
backward regions since the selectron exchange diagrams 
are dominant and give the (1 + cos^) 2 or (1 — cos 9) 2 an- 
gular dependence. We note that the contributions from 
the A3 = A4 case, which could give sin 2 9 dependence, are 
negligible for heavy selectron masses as mentioned before. 

Let us now turn to the simulation for the single-photon 
signal with missing energy. The partial decay rate of the 
neutralino decay into a photon and a gravitino is given by 
(see also Appendix |B|) 



I 2 m~ 
1 x{ 



487rM P1 ra 2 /2 



(19) 



and r(xi — > jG) — 0.21 GeV for the bino-like neutralino 
with m s o = 300 GeV and m 3/2 = 10~ 13 GeV. An ir- 
reducible SM background for the signal of mono-photon 
plus missing energy comes from e + e~ — > In addition 
to the minimal cuts for the detection of photons 



(20) 



6 K. Mawatari, B. Oexl, Y. Takaesu: Associated production of light gravitinos in e + e and e 7 collisions 



[fb] 


(P e -,P e+ )= (0,0) 


(0.9, 0) 


(0.9,-0.6) 




= 500 GeV 


mi = 400 GeV 


15 


23 


37 






800 GeV 


48 


75 


119 






1200 GeV 


64 


100 


159 






SM background 


1592 


178 


94 




= 1 TeV 


mg = 400 GeV 
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1443 
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Table 2. Cross sections in fb unit for the signal, e + e — > XiG — > 7GG, and the SM background, e + e — > 'yfu, at ^/s = 500 GeV 
and 1 TeV with different beam polarizations P e ±. We take m 3/2 = 10~ 13 GeV, m^o = 300 GeV, and B(x1 -> 7<5) = 1. The 
minimal cuts in (|20p and the Z-peak cut in (|21|l are taken into account. 



we impose the Z-peak cut 

™2 



0.02 



£ 7 < 



5r z , 



(21) 



which can remove the contributions from e + e~ — > 7Z — > 
71/^. The most significant background coming from the 
i-channel W^-exchange process can be reduced by using 
polarized e beams. 

In Table [21 the selection efficiencies for the signal and 
background processes with different polarization^ are pre- 
sented, where the above kinematical cuts, (|2T)]) and (|2"Tj) . 
are taken into account. The cross sections both for the sig- 
nal and background are calculated by MG/MEv4 [23] sup- 
porting gravitino interactions [TTll2"2"] . Here, we assume the 
branching ratio of the neutralino decay is unity, B(x± ~> 
7G) = 1, although other decay modes can be significant 
in some parameter space [21] . Since the cross section with 
e beam polarizations P e ± (\P e ± \ < 1) is given by 

. ( P e -,P e+) = 2£(I±f^)(I^). A , (22) 
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Fig. 5. Normalized energy distributions of the photon for 
e+e" -> x°G -> 7GG at ■s/s = 500 GeV (left) and 1 TeV 
(right), where rn s± = 400 (solid), 800 (dashed) and 1200 (dot- 
ted) GeV with m^o == 300 GeV are considered. The kine- 
matical cuts in (|20|) and (|21|) and the beam polarizations 
{P e - , P e +) = (0.9, —0.6) are taken into account. Those of the 
SM background are also shown by dot-dashed lines. 



the signal cross section for the bino-like neutralino can 
be enhanced, while the background can be reduced quite 
effectively, by using a positively polarized e~ beam (P e - > 
0) and a negatively polarized e + beam (P e + < 0). It must 
be noted again that the signal cross section is inversely 
proportional to the gravitino mass squared. 

Figure [5] shows normalized energy distributions of the 
photon for the signal and the SM background, correspond- 
ing to 20,000 events each, at y/s = 500 GeV (left) and 
1 TeV (right), where the selectron mass of 400, 800 and 
1200 GeV with the 300 GeV neutralino mass are consid- 
ered. The kinematical cuts in ([2"U| and (|2"Tj) and the beam 
polarizations (P e -,P e +) — (0.9,-0.6) are taken into ac- 
count. The signal distributions are flat, independent of the 
selectron massQ and restricted as 



■;?)- 



2VS <E " < 



(23) 



6 \Pe-\ > 0-8 and \P e +\ > 0.5 are designed at the Interna- 
tional Linear Collider (ILC) [23]. 

7 We point out that the photonic energy distributions of 
Fig. 5 in Ref. [8] should be flat and not depend on the se- 
lectron mass. 



where the lower edge can determine the mass of the neu- 
tralino. It should be stressed again that the photon dis- 
tribution is simply given by kinematical effects of the de- 
caying neutralino. The background is mostly distributed 
in the low-energy region, and hence a further cut on the 
photon energy would be useful to enhance the signal over 
background. 



Finally, we present the angular dependence of the pho- 
ton in the laboratory frame in Fig. [5] The original angular 
distributions of the neutralino in Fig. H] are flattened for 
the case of a 500 GeV collider since the neutralino decays 
isotropically in its rest frame and the boost effect is small. 
On the other hand, the angular distributions still survive 
for the case of a 1 TeV collider. This indicates a possi- 
bility to examine the mass of the t, it-channel selectrons 
when the decaying neutralino has a large momentum. We 
note that a kinematical cut on the forward and backward 
regions would also help to reduce the background. 
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Fig. 6. Normalized angular distributions of the photon in the 
laboratory frame for e + e~ — > XiG — > 7GG. The detail is the 
same as Fig. [5] 



Diagrams by MadGraph e- a -> er- gld 



graph 1 



graph 2 



N3 



. 3 



3 Selectron-gravitino production in e 7 
collisions 



graph 3 



graph 4 



In this section, we repeat a study as in Sect. [2] for a sce- 
nario of a selectron NLSP with a gravitino LSP. We con- 
sider associated gravitino productions with a selectron, 
especially a right-handed selectron, in e~7 collisions with 
the prompt selectron decay into an electron and a grav- 
itino, 



N4 



graph 5 




graph 6 



Fig. 7. Feynman diagrams for the e R -G production in e 7 
(24) collisions, generated by MadGraph 22 . 



leading to a mono-electron plus missing-energy signal. 



basis by 



3.1 Helicity amplitudes 



Here we present the helicity amplitudes explicitly for the 
production process: 



(pu y) + 7(Pa, A2) -> e£(p a ) + g( P4 , y ) . (25) 



\4 



The helicity amplitudes for the process are expressed as 
sums of s-, t-, and u-channel amplitudes 



iM 



AiA2,A4 



~ em f A 1 

= ~7^TT 1 e M<>2, A 2 ) 

V3MpiTO 3 /2 s 

x u{ Pil \ i )(f l +^ 2 )7 M u(pi,Ai), (27a) 

em^oC^'Clf i 

" 2V3 , Mpim 3/ 1 2 i~m|o ^ ^ 

x u(p4, A 4 )[|> 2 ,7 Al ](/ 1 -^ 3 +m x o)u(pi,Ai), 

(27b) 



-em- 



1 



V3Mpir71 3 /2 m - mg Ai 
x u(p 4 , A 4 )u(pi, Ai) (p 3 +pi - p 4 ) M , (27c) 



Ma 1 a 2 ,a 1 =M s +Y,M u +M u , 



(26) 



corresponding to the Feynman graph 6, (2+3+4+5), and 
1, respectively, in Fig. [7] They are given in the four-spinor 



with the couplings C 7 -^ and C± Xi defined in (fTU|) . Only 
the Ai = +1 case contributes to the (or e+ in our 
notation) production in the final state, which is relevant 
in our following analyses. On the other hand, the Ai = — 1 
case gives nonzero amplitudes only for the production. 
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Table 3. The reduced helicity amplitudes M\ 1 x 2 ,\ i for e Ai 7 Aa — > e R G Xi . 



We define the reduced helicity amplitudes, .A/f, as 



3.2 Cross sections and kinematical distributions 



iM Xl x 2M = j= C x/psMx lX2 Mi (28) 

\/6 A/pim 3/2 

and these are shown in Table [3] Similar to (TT2"1) . the four- 
momenta and helicities of the external particles are de- 
fined in the e^f CM frame with (3 = 1 — m§ /s. The fol- 
lowing features of the amplitudes are worth noting: 

1. The amplitude is zero since the coupling struc- 
tures do not allow this helicity combination for a mass- 
less goldstino. 

2. The overall angular dependence is dictated by J = 1/2 
or J = 1 d functions as 



M 



oc d 



IA1/2-A2I 

Ai/2-A 2 ,-A 4 /2 



(0). 



(29) 



3. The amplitude Ai* depends on the mass of the prop- 
agating neutralinos; as the neutralino mass increases, 

A4\. 1 _ becomes larger, while .M* , becomes smaller. 

On the other hand, the A4 S and Ai u do not depend 
on their mass but on the selectron mass. 

4. The right-handed selectron can couple only to the bino 
component of neutralinos i.e. Uu in the Xi~ e ~^R 
coupling in (fTU|) . Therefore, e.g. for the bino-like light- 
est neutralino case, only the Xi-exchange amplitude is 
nonzero among the four neutralino amplitudes. 

5. In the threshold region, similar to the e + e~ — > xjG 
process, the amplitudes are proportional to /3 3 / 2 , which 
gives rise to the strong suppression on the production 
cross section. 

6. In the high-energy limit, the amplitude .+ be- 
comes dominant, and hence the cross section depends 
on the neutralino mass but not on the produced selec- 
tron mass. 

We note that our helicity amplitudes in Table [3] agree 
with Eqs. (4) and (5) in [9]|f| where the heavy neutralino 
mass limit is assumed. The helicity-summed amplitude 
squared also agrees with Eq. (7) of [H] for the squark 
associated process qg —> qG after substitutions for the 
masses and the couplings as in (flUl) and the exchange of 
t «-» u. Moreover, we checked our amplitudes for each he- 
licity combination numerically by the gravitino/goldstino 
code in MG/ME [TTll2"2"]. 



8 Except the sign in the parentheses of the first term in 
Eq. (4) in [9]. 



In practice, a high-energy photon beam is provided by 
the backward Compton scattering of laser photons on a 
high-energy electron beam [24}[25] , as an option of a future 
linear collider [23] . Let us introduce the photon luminosity 
function [2"4ll2"g] 



fi(x,y) 



1 



N(y) 



1 



+ 1 



4r(l - r) 



1 - x 

+ P e P l ry(l-2r)(2-x) , (30) 



where x — E 1 /E e is the energy ratio of the scattered 
photon and the electron beam, y is a parameter controlled 
by the laser energy^ r — x/(l — x)y, and P e (Pi) is the 
incident electron beam (laser photon) polarization. The 
integral J f 1 (x, y) dx is normalized to unity by 



N(y)=(l Wl + y)+ + 

\ y y A J 2 

5 



+ PePl 



(l + -)ln(l + »)-2 + T 



1 



y 2(1 + y) 2 
1 



+ V 2(1 + yf 



(31) 



Figure [8] shows the luminosity function for different elec- 
tron beam polarizations, with y = 4.8 and Pi = —1.0 as an 
optimal parameter choice [28] . The maximal energy frac- 
tion is fixed by :r max = y/(l + y), e.g. x max ~ 0.83 for 



9 y = AE e Ei/m\ in the zero angle limit of the Compton 
scattering. 
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Fig. 8. The distribution functions of Compton back-scattered 
photons in (|30|) for different electron beam polarizations. 
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Fig. 9. Total cross sections of associated gravitino productions 
with a right-handed selectron in e~7 collisions, e~7 — > e^G, 
at ^fsZ = 500 GeV and 1 TeV for m 3/2 = 10~ 13 GeV as a 
function of the selectron mass. The neutralino mass is fixed at 
400 (solid), 800 (dashed) and 1200 (dotted) GeV, respectively. 
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Fig. 10. Angular distributions of the selectron in e 7 



at y/s = 455 (left) and 910 (right) GeV in the e-y CM frame, 
with 7713/2 = 10~ 13 GeV and = 300 GeV. The neutralino 
mass is taken to be 400 (top), 800 (middle) and 1200 (bot- 
tom) GeV, respectively. The photon-helicity dependent cross 
sections da^ are also shown by dashed lines (A2 = +1) and 
dotted lines (A2 = — 1). 



y = 4.8. The distribution with highly polarized electron 
and laser beams (P e Pi ~ —1) has a strong peak at the 
high-energy endpointF"! 

The full cross section at an e-f collider is calculated 
by convoluting the ej cross section (cr e7 ) with the photon 
distribution function of Eq. (f5U|) as 

<j(see) = f "" X Mx, y) a^(s) dx (32) 

with a; m i n = m\/s ee and s = xs ee , where ^/s ee is the 
original e~e~ CM energy. The spin-summed and averaged 
e7 cross section is obtained by a ei — {aV + erf_ 7 )/2 with 
the photon-helicity (A2) dependent cross section 

d °Z = Ysk^ ]MMX ^ m *- (33) 

Figure |H] shows total cross sections of the associated 
gravitino productions with a right-handed selectron in e~7 
collisions as a function of the selectron mass, where the 
CM energy, ^/s ee , of the e~e~ system is fixed at 500 GeV 
and 1 TeV. The parameters for the photon luminosity 
function in ((30)) are taken to be y = 4.8 and P e Pi = —0.9. 
For simplicity, we assume a bino-like lightest neutralino so 
that only the ^-exchange amplitude is taken into account 

10 Although the scattered photons are polarized when P e 7^ 
or Pi 7^ 25,28], we average the two opposite polarized 
modes, (P e > 0, Pi < 0) and (P e < 0,Pi > 0), so that we 
consider f-y(x,y) as the unpolarized distribution function in 
the following analyses. 



for the t-channel amplitude; see comment 4 in Sect. 13.11 
We note that a oc l/rriy 2 and we take m.3/2 = 10" 13 GeV 
in our study. When the selectron mass is close to the col- 
lider energy, the cross sections are strongly suppressed due 
to a ei oc /3 4 as mentioned in comment 5 in Sect. 13.11 
similar to the e + e~ — > process in Fig. [3] In addi- 
tion, the production cross section is nonzero only when 
Wen < V^maxSee- It should be stressed here that the cross 
section is quite sensitive to the mass of the ^-channel in- 
termediate neutralinos, even if the collider energy cannot 
reach their mass. 

Before we consider the selectron decay, let us look in 
detail at the angular distribution of the produced selec- 
tron in the ej system since the scalar decay — > eG is 
isotropic and hence the electron distribution is given by 
purely kinematical effects of the decaying selectron and 
the boost from the e~7 CM frame to the e~e~ laboratory 
frame. In Fig.[T0l the cos distributions of the selectron in 
e"7 e^G for m 3/2 = 10~ 13 GeV and m e - s = 300 GeV 
are shown, where the scattering angle 9 is defined from 
the momentum direction of the incident electron in the e"f 
CM frame. We fix the CM energy as y/s = ^x max s ee , 
where the photon luminosity sharply peaks for P e Pi ~ —1 
(see Fig.©, i.e. 455 GeV for a 500 GeV collider (left) and 
910 GeV for a 1 TeV collider (right). One can find that 
not only the total cross section as shown in Fig.[5]but also 
the angular distribution is quite sensitive to the mass of 
the i-channel intermediate neutralinos. As the neutralino 
mass is increasing, the cross section with A2 = +1 be- 
comes larger and the peak is shifted to the backward 
since the f-channel amplitude becomes more important 
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Cross sections in fb unit for the signal, e~7 — > 
e~GG, and the SM background, e~7 — > e~vi>, at 
= 500 GeV and 1 TeV without and with the electron 
beam polarization P e _ = 0.9. We take m 3 / 2 = 10 -13 GeV, 
rof H = 300 GeV, and B(e.R — > eG) = 1. The minimal cuts in 
(|35p and the Z-peak cut in (|3fc>p are taken into account. 



and its intrinsic sin | angular dependence is revealed as 
l/(t — to?) goes to 1/to?. On the other hand, the cross 
section with A2 = — 1 becomes smaller. As easily seen in 
Table O the productions to the forward region (cos# = 1) 
are forbidden for all helicity combinations because of the 
angular momentum conservation, while the productions 
to the backward region (cos 8 = — 1) are allowed only for 
the A2 = +1 case. 

Let us now turn to the simulations for the signal of 
single-electron plus missing energy in the e - e~ laboratory 
frame. The partial decay rate of the selectron decay into 
an electron and a gravitino is given by 



r{e R -> eG) 



71)'' 



48ttM p1 to 2 



(34) 



3/2 



and r(en - 
m 3 / 2 = 10" 



eG) = 0.27 GeV with m %R = 300 GeV and 
" 13 GeV. An irreducible SM background for the 
event of mono-electron plus missing energy comes from 
e~7 — > e~vv. In addition to the minimal cuts for the 
detection of electrons 



E e > 0.03^8, |?7 e |<2, 
we impose the Z-peak cut 



Mr, 



> 100 GeV, 



(35) 



(36) 



which can remove the contributions from e~7 — > e~ Z — > 
e~vv. The main background contribution coming from the 
W-exchange can be reduced by using a polarized electron 
beam. 

In Table |U the selection efficiencies for the signal and 
background processes without and with the electron beam 
polarization are presented, where the above two kinemat- 
ical cuts, (|3"5t and (|3l)|). are taken into account and it is 
assumed that the branching ratio of the selectron decay to 
an electron and a gravitino is unity, B(er — > eG) = 1. The 
cross sections both for the signal and background are cal- 
culated by MG/MEv4 [5T] with gravitino interactions [T71 
I22j . where we also implemented the photon luminosity 
function of ([3H| . By using a positively polarized electron 
beam of P e - — 0.9, the signal is enhanced by a factor of 
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Normalized energy distributions of the electron for 
e R G -> e~GG at ^7 = 500 GeV (left) and 1 TeV 
where mj = 400 (solid), 800 (dashed) and 1200 (dot- 
ted) GeV with mj a = 300 GeV are considered. The kinemat- 
ical cuts in ()35|) and ()36|) and the electron beam polarization 
P e - = 0.9 are taken into account. Those of the SM background 
are also shown by dot-dashed lines. 



1.9 because the cross section with Ai — — 1 is zero, while 
the background can be reduced by a factor of 10. It must 
be noted again that the signal cross section is inversely 
proportional to the gravitino mass squared. 

Figure [11] presents normalized energy distributions of 
the electron for the signal and the SM background, corre- 
sponding to 20,000 events each, at y/s^ = 500 GeV (left) 
and 1 TeV (right), where the lightest neutralino mass of 
400, 800 and 1200 GeV with the 300 GeV selectron mass 
are considered. The kinematical cuts in (1331) and (|36l) and 
the electron beam polarization P e - = 0.9 are taken into 
account. We notice again that the electron distribution is 
given by two boost effects, along the momentum of the 
decaying selectron and along the beam axis. The momen- 
tum of the incident electron is chosen to the +z direction, 
and hence the produced electrons in the e^ CM frame 
are boosted to the forward direction. Although the signal 
distributions no longer have either a flat shape or a sharp 
edge due to the boost along the beam direction, the energy 
is restricted as 



<E 7 < 



(37) 



where the lower edge can determine the selectron mass. 
The energetic electrons tend to be suppressed since the 
original selectron productions in the forward region are 
not allowed. Moreover, the z-axis boost effect makes the 
distributions slightly different for the different neutralino 
mass, which can be dictated by the peak shift in Fig. [TDJ 
Since the background is located mostly in the low-energy 
region, a certain amount of the energy cut can help to 
enhance the signal over the background. 

Finally, the angular distributions of the electron are 
shown in Fig. I12[ where the angle 9 e is measured from 
the direction of the electron beam, or the +z direction, 
in the e~e~ laboratory frame. The electrons tend to be 
produced more in the forward region (cos 8 e > 0) due to 
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Fig. 12. Normalized angular distributions of the electron in 
the e~ e~ laboratory frame for e~7 — > e^G — > e~GG. The 
detail is the same as Fig. 1111 



the system boost. For the ^fs^ = 500 GeV case the orig- 
inal dependence in the angular distributions of the 
selectron shown in Fig. [10] is no longer observed, while 
the dependence can be seen at y/s ee = 1 TeV. This in- 
dicates that we would be able to determine the mass of 
the ^-channel neutralinos when the collider energy is rel- 
atively higher than the selectron mass. We note that the 
background can be reduced further by a kinematical cut 
on the backward region. 

Before closing this section, we point out that selectron- 
neutralino associated productions in neutralino LSP sce- 
narios lead to the same signal, or e + and have been 
studied intensively [29,30,31,32^. Since the LSP mass is 
quite different between the two models, O(eV) for the G 
LSP and O(100 GeV) for the Xi LSP, the distributions of 
the final electron are distinctive and could provide a hint 
of SUSY breaking mechanism. 



4 Summary 

Associated gravitino productions with a SUSY particle 
can be observed in current and future collider experiments 
if the gravitino is very light. In this paper, we restud- 
ied the two associated-gravitino-production processes for 
a future e + e~/e~7 collider by using the gravitino imple- 
mented MadGraph/MadEvent. 

First, we studied gravitino productions in association 
with a neutralino which promptly decays into a photon 
and a gravitino at an e + e~ collider, e + e~ — > \iG — > 7GG. 
By using the effective goldstino interaction Lagrangian we 
explicitly presented the helicity amplitudes for the pro- 
duction process, which give us deep understanding for the 
threshold behavior and the SUSY-mass dependence of the 
production cross section and the angular distributions. We 
also examined selection efficiencies by kinematical cuts 
and beam polarizations for the signal and SM background 
processes, and showed that the energy and angular distri- 
butions of the photon in the final state can explore the 
mass of the i-channel exchange selectrons as well as the 
mass of the decaying neutralino. 



Second, we considered gravitino productions associ- 
ated with a selectron which subsequently decays into an 
electron and a gravitino at an ej collider, e~7 — > e~G — > 
e~GG. We repeated the same analyses as in the first pro- 
cess; we presented the explicit helicity amplitudes for the 
production process, and discussed the mono-electron plus 
missing-energy signal, including the energy spectrum of 
the backward-Compton scattered photons for incident pho- 
tons. Similar to the e + e~ — » x?G process, we found that 
the production cross section and the kinematical distribu- 
tions of the electron in the final state are quite sensitive 
to the mass of the i-channel intermediate neutralinos as 
well as the mass of the decaying selectron. 

We finally note that, throughout our study, we care- 
fully checked our calculations with the previous works 
both analytically and numerically, and pointed out a few 
disagreements. 

Before closing, we recall that all the helicity amplitudes 
we presented are easily applicable to qq — > gG and to 
qg — > qG subprocesses for hadron colliders. 
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A Effective goldstino interaction Lagrangian 

We briefly present the relevant terms of the interaction 
Lagrangian for our study. The effective goldstino interac- 
tion Lagrangian among goldstino, electron and selectron, 
G~e-e±, and among goldstino, neutralino and photon/Z- 
boson, G-Xi~V(= 7/^), in non-derivative form is 
9 

ITTt ~ 

V3Afpi m 3 / 2 

- Xi ^[ 7 ^,7"]^o(^K - d v V»\ (38) 

4V6Mpim 3/2 

where e± denotes the right-/left-handed selectron, P± = 
|(1 ± 75) is the chiral-projection operator, and the cou- 
pling C v ^ is defined in (TTUI); 

see more details in 1221 . All 

other relevant terms are 

Ceev = e^ e 7l4u " (9+ p + + 9-P-)Z^ e , (39) 



xe 
tively. 



(40) 
(41) 



where g± and G± Xi are defined in ^ and (fTU|) . respec- 
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B Helicity amplitudes for Xi — > lG 

We show helicity amplitudes for the neutralino decay into 
a photon and a gravitino, 



Xi(pi>y) -> 7(P2, A 2 ) + G?(p 3 , y 



A, 



(42) 



The partial decay rate in the neutralino rest frame is given 
by 



(43) 



-C* l m$o 

Al 



and the helicity amplitudes are calculated as 
M+,++ = -M 

M + 



cos — , 



v / 3Mpim 3/2 2 

-(77Xi m 3 o 

= M_, ++ = ^== ^sin-, (44) 



\/3Mpim 3 / 2 2 ' 

with C 7 * 1 in (fTU|) and the decay angle 0* defined from 
the quantization axis of the neutralino spin. The angular 
dependence is dictated by J = 1/2 d functions as 



1/2 

•MAi,A a Aa K d A 1 /2,A 2 -A 3 /2( 6 '*)- 



(45) 



Summing over the initial or final helicities for the ampli- 
tudes squared gives the isotropic decay distribution in the 
rest frame, and one can find the well-known decay rate in 
Eq. CE9]). 
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